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Abstract 

Self-similar properties of proton structure in the kinematic region of low Bjorken x are in- 
troduced and studied numerically. A description of the proton structure function F^x, Q 2 ) 
reflecting self-similarity is proposed with a few parameters which are fitted to recent HERA 
data. The specific parameterisation provides an excellent description of the data which 
cover a region of four momentum transfer squared, 0.045 < Q 2 < 120 GeV 2 , and of 
Bjorken x, 6.2 • 10" 7 < x < 0.01 . 



1 Introduction 



Recent measurements of the HI [[T|] and ZEUS [|2|] collaborations at HERA enable to study the 
proton structure in the region of low Q 2 < 1 GeV 2 where perturbative QCD has to face com- 
putation difficulties arising from the increase of the strong coupling constant a s (Q 2 ). Never- 
theless, there is a number of approaches to describe the transition to low Q 2 at small x together 
with the region, Q 2 > 1 GeV 2 , which is described by perturbative QCD very well. Such at- 
tempts involve Reggeon exchange ideas [|3|], dipole interactions [0j, vector meson dominance 
(VMD) [£p, efficient parametrisations [^] and others. 

This letter presents a different point of view, based on the idea that the proton structure at low 
x is of fractal nature. Using the fractal dimension concept [|7|], a simple parametrisation of the 
proton structure function F 2 (x, Q 2 ) is obtained with a few well defined parameters. A numeric 
study is made using recent small x HERA data, for Q 2 between 0.045 GeV 2 and 120 GeV 2 . 

2 Fractal Dimension 

The concept of fractal dimension requires to understand what is meant by dimension. In non- 
fractional dimensions a number of dimensions corresponds to a number of independent direc- 
tions in a corresponding coordinate system. For example, a line has obviously one dimension, a 
square two and a cube has three dimensions. The dimension of the Sierpinski gasket [^p, shown 
in Fig. HJ, needs a more general definition. 




Figure 1: Sierpinski gasket fractal in iterations No. 1, 3 and 6 (from left). Iteration No. 1 
corresponds to the seed image which is arbitrary while the iteration always converges to the 
same object. 

The cube, square and line are self-similar objects: when a line is broken in the middle two 
lines are obtained, each of half length. By magnifying one of them by a factor of two the original 
line is rebuilt. The same may be done by dividing a square into four small squares or a cube 
into smaller cubes. For example, when the magnification factor for a square is 3 the number 
of smaller squares will be 3 2 = 9, for a cube the same magnification gives 3 3 = 27 cubes. In 
general, when M is the magnification factor, then the number of objects will be M v , where V 
is the dimension of the object. The dimension V can thus be defined as 



V 



log M 
logM 



log(number of self-similar objects) 



(1) 



log(magnification factor) 
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According to this formula, the dimension of the Sierpinski gasket is fractional. When the 
magnification factor M is 2 there are 3 identical pieces of the gasket, for M = 4 the gasket 
consists of 9 small copies of itself. Therefore its fractal dimension is defined by 



v = log3 = log9 = l 5g496 
log 2 log 4 

Roughly speaking, the fractal dimension describes how complicated or how large a self- 
similar object is. A plane is 'larger' than a line. The Sierpinski gasket is not a line but also 
far from being a plane. Actually, there exist fractals which are constructed from lines but have 
dimensions 2 or 3, and therefore fill a plane or a space. An example of such a curve is the so 
called Hilbert curve. 

The definition of a dimension, given in equation (|l[), may be generalised for the case of 
non-discrete fractals. In this generalisation, the magnification (scaling) factor is a real number 
z and the number of self-similar objects is represented by a density function f(z). Taking into 
account that the dimension may change with scaling, a local dimension is defined as 



ologz 

For ideal mathematical fractals, discussed so far, V(z) is constant for the whole fractal. 
Introducing a scale dependent dimension is natural because many fractals in nature (e.g. plants 
or coastlines) are not mathematically ideal and usually have a fractal structure only for a certain 
region of magnification. In such a region, the dimension is approximately constant, T>{z) = V, 
and, following eq. (|3|), the density function f(z) is given as 



log/(» =V-logz + V (4) 

where T> defines the normalisation of f(z), which thus has a power law behaviour, f(z) oc z D . 

In general, fractals may have two independent magnification factors, z and y. In this case 
the density f(z,y) is written in the following way 



log/0, 2/) = V zy -log z - logy + V Z - log z + V y - logy + V . (5) 

Here the dimension V zy represents the dimensional correlation relating the z and y factors. 
The function f(z, y) satisfies a power law behaviour in z for fixed y and in y for fixed z. 

It is important to mention that there is certain freedom in selecting magnification factors 
without changing a shape of the function f(z, y). It is possible to use any non-zero power of a 
factor multiplied by a constant: z — > az x . The only effect of such a change is a redefinition of 
the dimensional parameters T^{ z ,y,zy} an d of the normalisation T> , respectively. 
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3 Self-similar Structure of the Proton 



Following the dimensional description of the presented fractal structures, it is interesting to 
study the properties of functions describing proton structure. In quantum chromodynamics the 
behaviour of the sea quark densities is driven by gluon emissions and splittings. The deeper 
the proton structure is probed, the more gluon-gluon interactions can be observed. These, in 
analogy to fractals, may follow self- similarity, i.e. scaling described by a power law. Indeed, 
there is a number of hints for a self-similar structure. As an example, Figure ^ shows the 
unintegrated u-quark density for fixed momentum transfer Q 2 and fixed Bjorken x, respectively. 




log(x) log(Q 2 ) 

Figure 2: Logarithm of the unintegrated u-quark density du(x,Q 2 )/dQ 2 as a function of 
Bjorken x (a) and Q 2 (b). The full and dashed lines correspond to GRV parametrisations in 
LO and NLO [§], respectively. 

For x < 0.01 (below the valence quark region) the unintegrated density function in log-log 
scale is linear. A linear behaviour is also exhibitted by the unintegrated density as a function of 
Q 2 for fixed x. Refering to eq. (Q), this suggests that x and Q 2 could be treated as appropriate 
magnification (scaling) factors. This is supporting the idea that the proton structure exhibits 
self-similar properties and may be described as a fractal object. 

Magnification factors are supposed to fulfil some criteria. They should be positive, non- 
zero and have no physical dimension. The two latter requirements concern the selection of Q 2 
as a magnification factor. The physical dimensionality may be removed by dividing Q 2 by a 
constant Qfo. For the case of Q 2 = 0, the non-zero requirement is not fulfilled, however, the 
access to this region is needed for integration of unintegrated densities. Thus instead of Q 2 a 
choice of 1 + Q 2 jQ\ as a magnification factor is appropriate. According to the freedom in the 
magnification factor selection, mentioned above, other equivalent choices are also possible, e.g. 
Ql/iQo + Q 2 )' (Qi + Q 2 )/lGeV 2 or similar combinations. It is also more appropriate to use 
1/x as a magnification factor rather than x itself: when the structure is probed deeper, x goes to 
zero while a magnification factor should rise. 
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4 Structure Function Parametrisation 



The concept of self-similarity, when applied to proton confinement structure, leads to a simple 
parametrisation of quark densities within the proton in a straightforward way based on Eq. (0). 
Using magnification factors 1/x and 1 + Q 2 /Q 2 ) , an unintegrated quark density may be written 
in the following general form 



log fi(x, Q 2 ) = V 1 ■ log - ■ log(l + %) + V 2 • log - + P 3 ■ log(l + %) + VI (6) 



where i denotes a quark flavour. Conventional, integrated quark densities qi(x,Q 2 ) are 
defined as a sum over all contributions with quark virtualities smaller than that of the photon 
probe, Q 2 . Thus fi(x, Q 2 ) has to be integrated over Q 2 , 



q i (x,Q 2 )= / f t (x,q 2 )dq 2 . (7) 



o 

Solving equation (|7|), the following analytical parametrisation of a quark density is obtained 

1 + P 3 -2? x log a: V Qo J 

Notice that in this parametrisation only the normalisation parameter T> l Q depends on the 
quark flavour while the other parameters are flavour independent. This assumption means that 
all quarks are following the fractal structure, i.e. the dimensions D L and the magnification 
factors are common for all of them and they differ in normalisation only. 

The proton structure function F 2 (x,Q 2 ) is related directly to the quark densities F 2 = 
x J2i e l(li + Qi)- Thus the assumption about the flavour symmetry of Eq. (||) allows to ex- 
press F 2 (x, Q 2 ) directly in the form given on the r.h.s. of Eq. (|S|) with x~ V2 replaced by x~ V2+l 
and with a common normalisation factor e v ° : 



FM = f S> (1 + - 1\ . (9) 

1 + P 3 - PxlogX V Qo J 

5 Fit to the Data 

The five parameters T>i and Qi are determined using recent data from the HERA experiments 
HI [0 and ZEUS [@] in the range 1.5 < Q 2 < 120 GeV 2 (HI) and 0.045 < Q 2 < 0.65 GeV 2 
(ZEUS). Additionally a cut x < 0.01 has been applied to exclude the valence quark region. The 
fit parameters are given in Table [j] and the corresponding description of the F 2 (x, Q 2 ) data is 
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Figure 3: Virtual photon-proton cross-section oy p oc F 2 (W 2 : Q 2 )/Q 2 as a function of Q 2 in 
W 2 bins. HI (points) and ZEUS (triangles) measurements are shown along with the fit to 4 
parameters (full line) and to all 5 parameters (dashed line). 
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T> 2 


£>3 


Q 2 [GeV 2 ) 


2 

X 


X /ndf 


all fit 


0.339 

±0.145 


0.073 

±0.001 


1.013 

±0.01 


-1.287 

±0.01 


0.062 

±0.01 


136.6 


0.82 


T>2 fixed 


0.523 

±0.014 


0.074 

±0.001 


1 

const. 


-1.282 

±0.01 


0.051 

±0.002 


138.4 


0.82 



Table 1: Results of the fit. The first row corresponds to a fit to all parameters, in the second row 
parameter V 2 was fixed to 1. The number of F 2 data points is 172, total errors were used for the 
X 2 calculation. 
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Figure 4: Measurement of the structure function F 2 (x, Q 2 ) as a function of x in bins of Q 2 by 
the HI experiment. The curve represents the fit to 4 parameters, which is indistinguishable from 
the 5 parameter fit in this kinematic region. 



shown in Figures [3], |] and [5|. The % 2 was calculated with total errors, adding the statistical and 
systematical errors in quadrature. When the relative normalisation of the HI and ZEUS data, 
which cover different Q 2 regions, was fitted no change beyond 1% was imposed by the fits. 
Thus the normalisations of the data sets were left untouched. 

Refering to Figure [3| the ratio F 2 (W 2 , Q 2 )/Q 2 is proportional to the virtual photon-proton 
cross-section a^* p (W 2 , Q 2 ). In the limit Q 2 — > and fixed W 2 the parametrisation @ behaves 
like Q 2 only for V 2 = I. This may be easily shown e.g. when the unintegrated structure 
function f(x, q 2 ) is introduced 
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rQ 2 

F 2 (x,Q 2 ) = / f(x,q 2 )dq 2 
Jo 



(10) 



the parametrisation of which is identical to @, with T> 2 replaced by T> 2 — 1 and V l replaced by 
V . If F 2 (x, Q 2 ) behaves like Q 2 for Q 2 — > then f(x, q 2 ) has to behave like a constant for any 
x = Q 2 1 (W 2 — M 2 ) — > 0. That is possible only if the divergent term, involving V 2 , is zero, i.e. 
for V 2 = 1. In this case, since other logarithmic terms go to zero, the ratio F 2 (W 2 , Q 2 )/Q 2 for 
Q 2 — ■> approaches the value eP° . 

In the fit with V 2 as a free parameter a value very close to 1 is obtained. Thus a second fit 
was made, where V 2 is fixed to 1 (see Table [j], second row). This fit has 4 parameters and gives 
nearly the same x 2 /ndf as the first fit to all 5 parameters. Within the kinematic range of the F 2 
data, both fits are nearly indistinguishable. As was stated above, the parameter V determines 
the virtual photon-proton cross-section in the photoproduction limit. Its value, obtained from 
the fit, gives 



a 



4.7r 2 a 



F 2 (W 2 ,Q 2 



189 ± 3 fib. 



(11) 



This is in approximate agreement with the total photoproduction cross-sections measured by 
the HI [p~0t] and ZEUS [O] collaborations which were not used in the fit. 
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Figure 5: Measurement of the structure function F 2 (x, Q 2 ) as a function of x in bins of Q 2 by 
the ZEUS experiment. The curve represents the fit to 4 parameters, which is indistinguishable 
from the 5 parameter fit in this kinematic region. 
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6 Summary 



The concept of the self-similar structure of the proton was introduced. This leads to a parametri- 
sation of the proton structure function F 2 (x, Q 2 ) which describes very well the low x HERA 
data, both in the non-perturbative and the deep inelastic domain. The introduced formalism 
uniquely defines the x and Q 2 dependence of parton densities, thus this approach is applica- 
ble also to other measures of proton structure, like the longitudinal structure function F L , the 
diffractive structure function F 2 D or the spin structure function gi. 
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